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ABSTRACT 


A three-dimensional , 54  degree  of  freedom  computer  model  of  the  heed/neck 
system  is  presented.  The  model  consists  of  a series  of  rigid  bodies  repre- 
senting the  bones  and  vertebrae  together  with  springs  and  dampers  representing 
the  muscles,  discs,  ligaments,  and  joints.  Equations  of  motion  sre  written  for 
the  model  by  using  Lagrange's  form  of  d'Alealbert's  principle.  Computer 
algorithms  are  developed  to  mmerlcally  calculate  the  coefficients  of  these 
equations.  The  governing  equations  are  then  Integrated  numerically  for  a 
number  of  specific  cases.  The  results  agree  very  well  with  experimental  data^ 
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NOMENCLATURE 

Disc  area 

Joint  damping  constants 
Disc  constants 
Standard  permutation  symbol 
Elastic  constant 

Disc,  ligament,  muscle,  or  Joint  force 
Generalized  force  associated  with 
Sheer  modulus 
Disc  height 

Second  cehtroldal  moment  of  area 

Centroldal  Inertia  matrix  of  the  kth  body  of  the  system 
Polar  moment  of  area 
Shape  factor 
Ligament  constants 

Muscle  constants  I J 

Mass  of  the  kth  body  of  the  system  Phi-]- 

> Unit  of  Vectors  fixed  in  R [FlllllLC 

■ Inertial  reference  frame  j f s 
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Partial  rate  of  change  of  position 

Generalized  coordinates 

Forward  direction 

Left  direction 

Axial  (up)  direction 

Displacement 

Rotation 

Partial  rate  of  change  of  orientation 


INTRODUCTION 

The  recent  research  efforts  in  biomechanics  and  biomedical  engineering  have 
fostered  a growing  interest  in  head/neck  biodynamic  models.  This  interest  has 
been  stimulated  by  the  awareness  that  60-70%  of  vehicle  related  fatalities  are 
caused  or  involve  injuries  to  the  head  or  neck.  To  date,  a satisfactory, 
comprehensive  head-neck  model  has  not  been  developed  which  accounts  for  all 
the  degrees  of  freedom  and  which  possesses  all  the  neuro-muscular  responses. 
There  have,  however,  been  a number  of  excellent  models  discussed  in  the 
literature.  Specifically,  in  1971,  Ome  and  Liu  (1J  developed  a discrete- 
parameter  spine  model  which  simultaneously  accounts  [or  axial,  shear,  and 
bending  deformation  of  the  discs,  for  the  variable  size  and  mass  of  the 
certebrae  and  discs,  and  for  the  natural  curvature  of  the  spine.  They  also 
present  an  extensive  literature  review  of  the  spine  models  prior  to  1970. 

Later  in  1971,  McKenzie  and  Williams  [2]  used  the  Orne-Liu  model  to  develop  a 
two-dimensional  discrete  parameter  head-neck-torso  model  for  "whiplash" 
investigation.  A two-dimensional  mechanical  linkage  model  simulating  head- 
neck  response  to  frontal  Impact  has  been  presented  by  Becker  [3).  This  model 
allows  for  elongation  of  the  neck.  It  concentrates  the  mass  at  the  head  mass 
center.  Springs  and  dampers  are  used  to  control  the  elongation  of  the  model. 

A three-dimensional  neck-torso  linkage  vehicle-occupant  model  has  been 
developed  by  Bowman  and  Robbins  [4].  The  model  has  two  ball-and-socket 
joints  and  the  neck  can  elongate  with  the  motion  limited  by  joint  stopping 
moments . 

In  addition  to  these  computer  models  there  have  also  been  developed  a 
number  of  wholebody,  anthropometric  dummy  models.  (These  are  currently  used 
extensively  by  the  automotive  industry).  In  1972,  Melvin,  et  al.  (5) 
presented  a mechanical  neck  for  anthropometric  dummies.  The  neck  consists  of 
three  steel  universal  joints  pinned  Into  aluminum  discs  with  shaped  rubber 
discs  around  the  joints.  The  joints  allow  the  neck  to  move  in  flexion, 
extension,  and  lateral  flexion  but  do  not  allow  for  either  rotation  or 
elongation.  A mechanical  neck  has  also  been  presented  by  Culver,  et  al.  [6]. 
It  consists  of  four  ball-joint  segments  and  one  pin-connected  "nodding" 
segment.  Viscoelastic  resistive  elements  Inserted  between  the  segments 
provide  for  bending  resistance  and  energy  dissipation  with  the  primary 
objective  being  to  model  flexion  and  extension  responses. 

In  this  paper,  there  is  presented  a comprehensive,  three-dimensional 
head-neck  computer  model  which  has  54  degrees  of  freedom  and  Includes  the 
effects  of  the  discs,  muscles  and  ligaments.  The  model  Is  based  primarily 
on  the  research  of  J.  Huston  and  Advani  [7,  8,  9)  which  in  turn  uses  the 
principles  of  the  whole-body  simulation  of  R.  Huston,  Passerello,  et  al. 

(10,  11,  12] 

The  balance  of  the  paper  is  divided  into  four  parts  with  the  first 
part  providing  some  preliminary  background  on  the  head-neck  anatomy,  geometry, 
and  modelling.  The  next  part  presents  a summary  of  the  development  of  the 
equations  of  motion  and  the  computer  algorithms  of  the  head-neck  model.  The 
third  part  presents  a comparison  of  numerical  results,  predicted  by  the  model, 
with  experimental  data.  The  final  part  contains  a brief  discussion  and  some 
conclusions . 
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HEAD-NECK  ANATOMY  AND  MODELLING 


A comprehensive  presentation  of  the  head-neck  anatomy  may  be  found  In 
references  [13-20],  The  anatomy  Is  conveniently  divided  Into  two  categories: 
bones  and  soft  tissue. 


Bones 

The  largest  and  heaviest  bone  Is  the  skull  which  consists  of  a large  cranial 
cavity  (enclosing  the  bran)  and  smaller  bones  containing  the  face  and  Jaw. 

The  skull  Is  actually  composed  of  21  closely  fitted  bones.  The  other  bones 
of  the  head-neck  system  are  seven  cervical  vertebrae  (C1-C7)  which  support 
and  provide  mobility  ot  the  head.  The  first  of  these  Cl,  called  the  "atlas", 
supports  the  skull.  The  second  C2,  called  the  "axis".  Is  distinctive  because 
of  Its  adontold  process  (or  axis)  which  rises  perpendicularly  to  the  vertebrae. 
The  five  remaining  cervical  vertebraw  are  roughly  annular  in  shape  and  are 
similar  to  each  other  with  a slight  Increase  In  size  going  down  from  C3  to  C7. 


Soft  Tissue 

The  soft  tissue  is  composed  primarily  of  the  discs,  the  muscles,  the  ligaments 
and  the  brain.  The  discs  provide  the  cushioning  or  separation  for  the 
vertebrae.  They  are  annular  in  shape.  The  ligaments  connect  the  cervical 
vertebrae  to  each  other  and  thus  allow  for  the  gross  and  fine  movement  of  the 
head  and  neck.  The  movements  of  the  head  and  neck,  which  may  be  classified 
grossly  as  flexion,  extension  and  rotation,  are  controlled  by  the  muscles. 

The  muscles  originate  on  the  various  cervical  vertebrae,  the  skull,  the  spine 
and  the  shoulder  bones.  The  brain  tissue  Is  basically  four  mass  volumes 
composed  of  two  cerebral  hemispheres  In  the  upper  half  of  the  skull,  the 
triangular  shaped  cerebellum  in  the  lower  posterior  and  the  brain  stem  In  the 
center  of  the  skull. 
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Modelling 

The  head-neck  system  is  modelled  by  a system  of  9 rigid  bodies  representing 
the  skull,  vertebrae  and  torso,  as  shown  in  Figure  1,  and  springs  and  dampers 
representing  the  discs,  ligaments  and  muscles.  The  masses,  inertia  matrices, 
and  overall  geometry  of  the  rigid  bodies  are  adjusted  to  match  the  actual 
human  values  [7].  Each  body  has  6 degrees  of  freedom  and  hence  the  entire 
system  has  a total  of  34  degrees  of  freedom. 

Following  Orne  and  Liu  [1],  the  discs  are  modelled  In  the  axial 
direction  as  two-parameter  viscoelastic  solids  with  the  uniaxial  force- 
displacement  relationship  being: 

F - (A/h) (dj«  + d2«)  (1) 


In  bending  and  shear  the  discs  are  modelled  as  linear  elastic  solids.  Using 
the  principles  of  strength  of  materials  theory  [7],  the  following  force  and 
moment  equations  are  developed: 
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where  as  shown  In  Plgure  1,  Z Is  In  the  axial  (up)  direction,  X Is  forward 
and  T la  to  the  left. 

The  ligaments  are  modelled  as  non-linear  elastic  bands  capable  of 
exerting  force  only  In  tension.  The  force-displacement  relation  Is  taken  as: 


F - IjS  + l2S 


(10) 


The  muscles  are  modelled  as  two  parameter,  visco-elastic  solids,  which 
like  the  ligaments,  only  exert  force  when  In  tension.  The  force-displacement 
relation  is  taken  as: 


F - m1«  + m2« 


(11) 


The  Joint  constraints  (limiting  the  relative  motion  of  the  bodies)  are 
modelled  as  one-way  dampers.  The  force-displacement  and  moment-rotation 
relations  are  taken  as: 


p ^ -C4  for  4X5 
0 for  4<0 


and  . . (12) 

-C8  for  8>0 

M 0 for  0<O 
where  the  damping  constant  Is 
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C + C. (X-Xmax)  for  X>Xmax 
C - C°  for  Xmln<X<Xmax 

C°  + C, (X  Xaln)  for  X<Xmln 

O 1 

where  X,  Xmax,  and  Xmln  are  the  values  of  the  displacement  or  rotations 
variable  and  Its  corresponding  maximum  and  minimum  values. 

The  values  of  these  various  constants  for  the  discs,  ligaments, 
muscles  and  Joints  for  the  various  directions  and  motion  are  difficult  to 
specify  precisely  due  to  a lack  of  experimental  data.  However,  the  values 
for  the  discs  may  be  obtained  from  Markolf  and  Steldal  [21],  Orne  and  Liu  [1], 
and  McKenxle  and  Williams  [2].  The  ligament  and  muscle  attach  points  may  be 
obtained  from  Francis  [22],  Lanier  [23],  and  Todd  and  Llndala  [24],  with  the 
spring  and  visco-elastic  constants  obtained  from  Nunley  [25]  and  Close  [26]. 


EQUATIONS  OF  MOTION 

The  model  of  Figure  1 can  be  considered  as  a "general  chain  system" — that  Is, 
a system  of  rigid  bodies  arbitrarily  assembled  such  that  no  closed  loops  are 
formed.  (Such  systems  are  sometimes  called  "open-chains"  or  "open-trees".) 
Using  Lagranges  form  of  d'Alembert's  Principle  as  developed  by  Kane  et  al. 
[27-30],  R.  Huston  and  Fasserello  [31,32]  have  developed  procedures  and 
computer  algorithms  to  study  such  systems.  Indeed,  they  have  shown  that  the 
equations  of  motion  may  be  written  In  the  form: 

ai]*J  " fl  (1,  J - 1 H)  (1*) 

where  N Is  the  number  of  degrees  of  freedom  of  the  system,  X,  are  the 
generalized  coordinates,  and  a,,  and  f are.  In  general,  functions  of  X, 
and  X,.  (Regarding  notation,  unless  otherwise  specified,  there  Is  a sum 
over  repeated  Indices  over  the  range  of  that  Index). 

For  the  system  shown  in  Figure  1 the  number  of  degrees  of  freedom  Is 
54,  and  hence  54  generalized  coordinates  (27  translation  and  27  rotation) 
need  to  be  defined  to  describe  the  configuration  of  the  system.  It  Is 
convenient  to  define  these  coordinates  as  relative  or  local  coordinates  (l.e. 
between  the  respective  bodies)  as  opposed  to  absolute  or  global  coordinates 
(l.e.  with  respect  to  Inertia  space).  Upon  Introducing  these  coordinates, 
the  angular  velocities  relative  to  an  Inertia  frame  of  the  various  bodies  of 
the  system  may  be  written  In  the  form: 

Sfk  - (1  - 1,  2,  3;  J - 1 27;  (15) 

1Jk  J"1  k ■ 1, . . . ,9) 

where  N.  are  mutually  perpendicular  unit  vectors  fixed  In  an  Inertial 
reference  frame  R and  where  the  i». are  called  "partial  rates  of  change  of 
orientation”  [29]  and  are  in  general  functions  of  X..  Similarly,  the 
velocities  of  the  mass  centers  In  R of  the  various  bodies  of  the  system  may 
be  written  In  Che  form: 

vk  - v^X^  W - 1.  2,  3;  J - 1 54;  (16) 
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where  the  v...  are  called  "partial  rates  of  change  of  position”  [29]  and  are 
In  general  functions  of  X..  Langranges  form  of  d'Alembert's  principle  then 
leads  to  governing  equations  of  the  form  of  Equations  (14)  with  N equal  to  54, 
and  with  a and  f ^ given  by  the  expressions: 


and 


*ij  mpVpJkVplk  + Ipknupjn“pik 

fl  ^F1  + BpVplkVpqkXq  + Ipkn“pln“pqkXq 


(17) 

(18) 


+ e,iu  u iii  ..I  XX) 
nmk  pqn  par  plk  pmr  q s 


where  F represents  the  generalized  force  with  respect  to  X and  it  may  be 
expressed  as  (27]:  1 


F.  * v ..F  . + u M 
1 plk  pk  plk  pk 


(19) 


where  F and  M ^ are  the  N.  components  the  equivalent  centroldal  force  and 
torque  or  the  equivalent  couple  acting  on  the  pth  body  of  the  system. 

Equations  (14)  then  represent  a system  of  54  simultaneous,  second 
order,  ordinary  differential  equations.  They  may  be  Integrated  numerically 
using  a standard  Integration  scheme.  Results  of  such  an  Integration  for  a 
variety  of  specific  cases  are  presented  In  the  next  part  of  the  paper. 


EXPERIMENTAL  VERIFICATION 


Several  experiments  have  been  conducted  which  may  be  used  to  obtain  a 
validation  of  the  model.  In  one  of  these,  a seated  cadaver  was  subjected  to 
head  Impacts  by  a rigid  pendulum.  Accelerometers  were  used  to  measure  the 
resultant  frontal  and  occipital  head  Impact  forces  and  accelerations.  Using 
the  Impact  force  data  as  Input,  the  acceleration  was  calculated  using  the 
computer  model.  A comparison  of  the  results  for  two  of  the  frontal  Impact 
experiments,  6-2  and  6-5  la  shown  in  Figures  2-5. 

In  the  same  set  of  experiments,  high-speed  cameras  were  used  to 
measure  the  acceleration,  velocity,  and  displacement  of  the  mass  center.  A 
comparison  of  the  results  with  those  predicted  by  the  computer  model  for 
experiments  6-1  and  6-2  are  shown  In  Figure  6. 

Finally,  the  model  was  checked  against  live  human  data  generated  by 
Ewing  and  Thomas  [33]  using  elaborate  testing  facilities.  A comparison  of 
the  results  for  the  head  angular  acceleration,  angular  velocity,  and  angular 
displacement  Is  shown  In  Figures  7,  8 and  9. 
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Fig.  6 Comparison  of  model  and  experiment  for 
head  mass-center  displacement,  velocity,  and  acceleration 


Fig.  9 Comparison  of  model  and  Ewing 
experiment  for  head  angular  displacement 


CONCLUSIONS 

Figs.  1-9  show  that  there  is  a very  good  comparison  between  the 
experimental  results  and  those  predicted  by  the  computer  model.  This 
is  Indeed  encouraging  and  it  suggests  that  the  model  represents  one  of 
the  most  sophisticated  head-neck  models  available.  However,  more 
testing  and  refining  needs  to  be  done.  Specifically,  the  three- 
dimensional  features  of  the  model  need  to  be  further  checked  with 
experimental  data.  Also,  better  experimental  values  for  the  soft 
tissue  mechanical  properties  needs  to  be  obtained.  Finally,  better 
criteria  for  injury  need  to  be  established.  When  this  is  done,  the 
model  will  provide  an  even  more  effective  and  economical  tool  for 
predicting  injury  in  a variety  of  high-acceleration/accident-configura- 
tion environments. 
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